SINKING OF ^^NE IN LIQUID WD INTERIORS 
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, Abstract. We assess the impact of the trace element ^^Ne on the coohng 

On \ and seismology of a liquid C/0 white dwarf (WD). Due to this elements' 

' neutron excess, it sinks towards the interior as the liquid WD cools. The 

\ gravitational energy released slows the WD's cooling by 0.5-1.6 Gyr. In 

^ ' addition, the ^^Ne abundance gradient changes the periods of the high 

p . radial order g-modes at the 1% level. 
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1. Biasing of ^^Ne Diffusion in Liquid White Dwarf Interiors 

After -"^^C and ""^^O , the most abundant nucleus in a M < Mq white dwarf 
■ (WD) interior is ^■^Ne . As the slowest step in the CNO cycle is the p+^^N 

, reaction, almost all CNO nuclei end up as ^^N at the end of H burning. Dur- 

ing the He burning stage, the reaction ^^N (a,7)^^F (/3+)^^0 (Q,7)^^Ne 
processes all ^^N into ^^Ne . For recently formed WDs of M < Mq this 
results in a ^^Ne mass fraction X22 ~ ^CNO ~ 0.02 (see, for example, [9]). 

^^Ne has a mass to charge ratio, A/Z greater than that of "^^C and 
^^O . As pointed out by [4] and [2], this affects the dynamics of ^^Ne nuclei 
in the following manner. In the degenerate WD interior, there exists an 
upward pointing electric field of magnitude eE k, 2mpg [2], g being the local 
gravitational acceleration. The resulting net force on ■^^Ne , F = —2mpgr, 
biases ^^Ne 's diffusion inward. In contrast, the predominant ^^C and ^^O 
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ions experience no net force. The resulting sedimentation of ^^Ne will 
impact the cooling history and seismology of the WD. 
The ^^Ne evolution is governed by mass continuity 

^ + V.J.. = 0, (1) 

where P22 and J22 are the local ^^Ne density and flux respectively. The 
flux is given by 

J22 = (-^^ - Vp22)i , (2) 

where D is ^^Ne's diffusion coefficient in the background plasma and v 
is the magnitude of its local drift velocity. D and v are related via the 
Stokes-Einstein equation: v = 2mpgD/kT. 

As discussed in [2] and [6], D is not well known for conditions in 
the WD interior. So, for this work, we took as a starting point the self 
diffusion coefficient of the one-component plasma, Dg, as calculated in 
[8],Ds ^ Sujpa^r-'^/^, where ioj = 4TTni{Ze)^ /ArUp, F = {Ze)'^/akT and 

= SAnip/Airp is the inter-ion spacing of the background ions. At F = 173, 
the plasma crystallizes [7]. We quantify the effects of the uncertainty of D 
on our results by performing calculations with varying D. 

2. ^^Ne Sedimentation and WD Cooling 

We simultaneously evolve the initially constant X22 and WD core temper- 
ature, Tc. The ^^Ne evolution is modeled as above. We halt sedimentation 
in the solid phase. We consider the WD interior to be isothermal, and 
calculate Tc{t) using 

—JT / -A dm = -L + Lg + , (3) 

dt Jo Amp ^ Amp dt 

where Cy is the heat capacity, L is the WD luminosity, Lg the power gener- 
ated by ^^Ne sinking, I is the latent heat released per nucleon, and Mcrys is 
the mass of the WD that has crystallized. L is determined from the L — Tc 
relations of [1] and Lg = F{J22/ P22)n224^TTr^dr, where 77-22 = P22 /22mp. 
We ignore fractionation effects. Finally, throughout we utilize background 
WD models composed of a single ion species, representing realistic C/0 
WD compositions by taking A = 14 and Z = 7. 

The overall impact of ^^Ne sedimentation on the WD age depends on 
the diffusion rate and when the WD crystallizes. Sinking is fast in massive 
WDs, but crystallization occurs earlier, somewhat mitigating the impact. 
The opposite is the case for low mass WDs, where sedimentation is slower, 
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Figure 1. Cooling of a O.6M0 WD. The lower dashed curve shows the results neglecting 
^^Ne . The remaining curves show (from bottom to top) the results for D = Q, Da, 5Ds, 
and lODs for the specified X22 value. The shaded regions gives the age range of NGC 
6791. 
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Figure 2. Theoretical CMD for NGC 6791, showing the cooling tracks for 0.6 and I.OM0 
WDs. Locations of WDs today for (bottom to top) D = 0, Ds, 5Ds, lODs are shown by 
the dots. 
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but can continue for up to 10 Gyrs. Ne content is also a factor, as Lg 
depends on the X22. Figure 1 shows cooHng curves for a O.6M0 WD model, 
comparing differing values of D and X22- The higher metallicity is that of 
NGC 6791, a metal rich, old (8.0±0.5 Gyr) open cluster [5]. Figure 2 details 
the expected CMD location of WDs in NGC 6791 today assuming several 
values of D, demonstrating how observations can constrain on the value of 
D. 

3. Enhanced Internal Buoyancy from ^^Ne 

The gradient in ^^Ne abundance creates a gradient in the mean molecu- 
lar weight per electron, /Xej that provides an extra restoring force for the 
g-modes in the WD interior. We can approximate the Brunt- Vaisala fre- 
quency, A^^ , in this case by considering an isothermal plasma of degenerate, 
relativistic electrons and ideal ion gas with a ji^ gradient due to X22 gra- 
dients. Doing so, we find (see [6] for details) 

h'^Hi'mp 11 dz 

where h = P/pg and ^Uj is the ion mean molecular weight. The first piece is 
the thermal contribution to N"^ while the second is the ^^Ne contribution. 

We estimate the impact of ^^Ne on the mode frequencies by performing 
the integral / Nh/rdilnP) through the liquid WD interior (the envelope's 
contribution to the integral is « 0.3 rad s~^). Neglecting the ^^Ne contri- 
bution to iV^ gives 0.073 rads~^ through the WD core, while taking into 
account the ^^Ne gradient gives 0.079(0.095) rad s"^ for D = Ds{5Ds) . 
Including the envelope contribution shows ^^Ne can impact mode periods 
by more than 1%, a level over 100 times larger than the measurement er- 
rors (see, e.g., [3]). Clearly, if meaningful results are to be inferred from 
WD pulsation studies, ^^Ne gradients need to be taken into account. 
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